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We report on a precise computation of the topological charge distribution in the SU(3) Yang-Mills 
theory. It is carried out on the lattice with high statistics Monte Carlo simulations by employing 
the definition of the topological charge suggested by Neuberger's fermions. We observe significant 
deviations from a Gaussian distribution. Our results disfavour the 6 behaviour of the vacuum energy 
predicted by instanton models, while they are compatible with the expectation from the large Nc 
expansion. 
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Introduction. — The 9 dependence of the vacuum en- 
ergy P, 0, d, 01 , or equivalently the functional form of the 
topological charge distribution, is a distinctive feature of 
the ensemble of gauge configurations that dominate the 
path integral of a Yang-Mills theory. In the Euclidean 
space-time the ground-state energy F{9) is defined as 



(1) 



where, as usual, (...) indicates the path-integral aver- 
age (our normalization is i^(0) = 0). In the large volume 
regime F{0) is proportional to the size V of the system, a 
direct consequence of the fact that the topological charge 
operator Q is the four-dimensional integral of a local 
density. The function F{9) is related to the probabil- 
ity of finding a gauge field configuration with topological 
charge Q — v hy the Fourier transform 



Pu = 



(2) 



Large Nc arguments [SJ, with Nc being the number of 
colors, suggest that the fluctuations of the topological 
charge are of quantum non-perturbative nature [y^ 0] ■ 
The 9 dependence of the vacuum energy is expected at 
leading order in l/Nc, and the normalized cumulants 
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which should scale asymptotically as Nc~^"' [a,LZll) have 
to be determined with a non-perturbative computation. 
On the other hand several models, such as the dilute gas 
or liquid of instantons, assume that the path integral is 
dominated by semiclassical configurations [1, H, [13, [HI ■ 
They predict a 9 behaviour of the form 



j^inst^^) = -VA{cos{9) - 1} , 



(4) 
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with A being exponentially suppressed at large Nc. 

The 9 dependence of the vacuum energy plays a crucial 
role also in the solution of the so-called U{1)a problem 
in QCD. The Witten-Veneziamo mechanism relates the 
cumulants of the topological charge distribution in the 
Yang-Mills theory with the leading anomalous contri- 
bution to the mass and scatte ring amplitudes of the 77' 
meson in QCD [H, [H, [11 [11, [ij. The known value of 
Ci in the SU(3) theory supports indeed the fact that the 
bulk of the rj' mass is due to the anomaly Tt'I . 

Recent theoretical developments in lattice gauge the- 
ory made it possible to find an unambiguous definition 
of the topological charge distribution with a finite and 
universal continuum limit 15, 18, 19, 20, 21, 22]. The 
aim of this work is a precise computation of the distri- 
bution of the topological charge in the SU(3) Yang-Mills 
theory. We observe significant deviations from a Gaus- 
sian behaviour: they disfavor the 9 dependence given in 
Eq. (HI) , while they are compatible with expectations from 
the large Nc expansion. 

In the past the distribution of the topological charge 
was already studied (see Ref. QJ^ and references therein). 
These computations, however, were not precise enough to 
observe deviations from the leading Gaussian behaviour. 
In this respect we have exploited the efficiency of the al- 
gorithm for the determination of the charge developed 
in Ref. [1^. Properties of the charge distribution have 
been investigated also with fermionic and bosonic meth- 
ods (see Refs. [H, [2^ and references therein). These 
results, however, are affected by model-dependent sys- 
tematic errors that are not quantifiable, and their inter- 
pretation rests on a weak theoretical ground. 

Topological charge definition. — The Neuberger-Dirac 
operator D is defined as 
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where I?w is the standard Wilson-Dirac operator and s 
is an adjustable parameter in the range |s| < 1 (for no- 
tations not explained here see Ref. [23|). It satisfies the 
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Ginsparg- Wilson relation ^2d\, and therefore the associ- 
ated fermion action preserves an exact chiral symmetry 
at finite lattice spacing [l^. The corresponding Jaco- 
bian is non-trivial, and the chiral anomaly is recovered d 
la Fujikawa [l^ [23] with the topological charge density 
operator defined as [23| 



(7) 



where the trace runs over spin and color indices. With 
this definition the topological charge in a given back- 
ground is given hy Q = J2xli^) ~ n+ — n^, with n+ 
(n_) being the number of zero modes of D with posi- 
tive (negative) chirality. The normalized cumulants C„ 
are thus defined as the integrated connected correlation 
functions of n charge densities^ 



= IT ^ {q{xi)...q{x2n)Y 

Xi,...,X2n 



(8) 



They have an unambiguous finite continuum limit which 
is independent of the details of the regularization [isl . 
[2ll . \2^ . At finite lattice spacing they are affected by 
discretization errors which start at O(a^). 

The large volume limit. — Being v integer- valued, F{9) 
is a periodic function with period 2tt. In the interval 
— TT < < TT it has its absolute minimum at 9 = 0, and 
it may then be expanded as 



n+1 



(2n) 



(9) 



Leading finite-size effects in the C„ are exponentially 
suppressed at asymptotically large volumes. They are 
proportional to e~^^g^, with Mg ~ 1.6 GeV being the 



lightest glueball mass [28| , and they become rapidly neg- 
ligible as soon as L is larger than 1 fm or so [29]. By 
inserting Eq. © in Eq. and neglecting exponentially 
small corrections proportional to e~ , we can express 
the topological charge distribution at large volumes 
by a saddle point expansion (usually named Edgeworth 
expansion in statistics^). This is an asymptotic series in 
powers of 1/V (or 1/{N^V) according to large Nc) which, 
up to higher order corrections, reads 



P. = 



V27 



l + ^^He4(;./a) 



(10) 



The parameters are — VCi and r = C-^JG ^ , a nd the 
Hermite polynomial He4 can be found in Ref. |31j . 



^ Correlation functions of an odd number of topological charges 
vanish thanks to the invariance of the theory under parity. 
The Edgeworth expansion is usually adopted in the context of 
the central limit theorem [30II . In our case the volume y, or 
(N'^V) at large A^c, plays the role of the number of independent 
degrees of freedom. 



The semiclassical models provide a sharp prediction of 
the topological charge distribution. By inserting Eq. ^ 
in Eq. ([2]), and taking into account that v is an integer, 
we obtain 



Inst 



e-^^/„(FA) 
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where are the modified Bcssel functions of the first 
kind [3l|. By construction all normalized cumulants are 
equal to A. 

Lattice computation. — The numerical computation is 
performed by standard Monte Carlo techniques. The en- 
sembles of gauge configurations are generated with the 
Wilson action and periodic boundary conditions. Each 
update cycle consists in 1 heat-bath and several over- 
relaxations of all link variables (more details can be 
found in Ref. [1^). The charge density is defined as in 
Eq. ([7]) with s = 0.4, and the corresponding topologi- 
cal charge has been computed by counting the number 
of zero modes of the Neuberger-Dirac operator with the 
algorithm proposed in Ref. [23| . 

The list of lattices, the value of the bare coupling con- 
stant P = 6/(7q, the linear size L/a in each direction, and 
the number of independent configurations are reported in 
Table m Lattice spacings and volumes have been chosen 
to have normalized cumulants with small discretization 
and finite-size errors. To estimate discretization effects 
we have simulated three lattices, A1-A3, with the same 
physical volume but different lattice spacings. Two addi- 
tional lattices, Bi and Ci, have been generated to quan- 
tify the magnitude of finite-size effects: they have the 
same bare coupling of Ai but larger volumes. The au- 
tocorrelation function of the topological charge has been 
computed at the three values of /3 by monitoring its value 
for several thousands of consecutive update cycles of the 
A lattices. The corresponding autocorrelation time is 
between 20 and 40 cycles. Since the cost of a Monte 
Carlo update is negligible with respect to the computa- 
tion of the index of the Neuberger-Dirac operator, we 
separated subsequent measurements on all lattices by a 
number of update cycles between 1 and 2 orders of mag- 
nitude larger than the above estimates. Statistical errors 
are then computed with the jackknife method by consid- 
ering the measurements as independent. A preliminary 
analysis of a subset of our results was presented at the 
conference "Lattice 2006" [32j . 

The Monte Carlo technique adopted here generates the 
gauge configurations with a probability density propor- 
tional to e~'^^", with Sym being the chosen discretiza- 
tion of the Yang-Mills action. This algorithm performs 
an importance sampling of the topological charge with 
the probability distribution given in Eq. A statis- 
tical signal for the n"^ cumulant is then obtained only 
if the number of configurations is high enough for the 
sample to be sensitive to terms suppressed as V^~^ in 
the Edgeworth expansion. For instance, the estimators 
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Lat 13 


L/a 


ro/a 


L[fm] 


N 


'2 IN orm 
Xdof 


2 inst 
Xdof 


iJ ,ii)clgc 
Xdof 


Ai 6.0 


12 


5.368 


1.12 


34800 


15 


27 


1.5 


A2 6.0938 


14 


6.263 


1.12 


30000 


12 


34 


1.4 


A3 6.2623 


18 


8.052 


1.12 


30000 


13 


41 


1.1 


Bi 6.0 


14 


5.368 


1.30 


30000 


1.3 


6.7 


0.14 


Ci 6.0 


16 


5.368 


1.49 


10000 


0.67 


2.4 


0.79 



TABLE I: Simulation parameters, number A'^ of configura- 
tions generated, and values of Xdoi for tlie fit of the data to 
a Gaussian (Norm), to the instanton prediction in Eq. (Illf) 
(Inst) and to the Edgeworth expansion in Eq. (|10p (Edge). 



of the first two cumulants 
1 ^ 
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with Vi being the value of the topological charge for a 
given gauge configuration and A'^ the total number of 
configurations, have variances which, up to sub-leading 
corrections, are (2(7" + a'^T)/N and (24cr^ + 72(t^t)/N 
respectively. The number of configurations on our main 
set of lattices, the A series, has been fixed to have a pre- 
cision of 15 — 20% on the second cumulant C2. For the 
lattice Bi the number of configurations is chosen so to 
have a signal for (Q'*)'^"" and therefore a rough estimate 
of finite-size effects^. Lattice Ci has been simulated to 
quantify finite-size effects in (Q^) with confidence. 

For each lattice we have compared the histogram 
of the topological charge distribution with three func- 
tional forms: a simple Gaussian centered at the origin, 
the Edgeworth expansion in Eq. (|10p , and the prediction 
from instanton models in Eq. (jlip . The free parameter(s) 
of each function has(ve) been fixed by maximizing the 
likelihood. For the symmetrized histograms the values of 
per degree of freedom at the minimum are reported 
in Table U and for the lattice A2 the data points and the 
three curves are shown in Fig. [1] The Edgeworth expan- 
sion reproduces well the behaviour of the numerical data 
at these volumes and lattice spacings within our statisti- 
cal errors*. On the A lattices the Gaussian distribution 
is incompatible with the data, while at the two larger vol- 
umes the Xjqj is still rather good. The functional form 



^ Given the scaling of the statistical error with V and n, it is very 
inefficient to compute higher cumulants at large volumes with the 
standard sampling procedure. Once is known, however, one 
could devise an adaptive importance sampling algorithm which 
integrates this information so to have a reduced o-^g. 
We also fitted the data with the functional form in Eq. (|2} and 
two non-vanishing cumulants. The conclusions are analogous. 
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FIG. 1: Number of configurations vs the topological charge 
for the lattice A2. To guide the eye, lines connect the values 
of the fitted distributions: the simple Gaussian (Norm), the 
instanton prediction in Eq. (|ll|l (Inst) and the Edgeworth 
expansion in Eq. (|10|l (Edge). The plots on right are a blowup 
of the top and the bottom of the distribution. 



suggested by instanton models is excluded on the lattices 
A1-A3 and Bi, and is off by more than two sigmas on 
the lattice Ci. On the A lattices a fit limited to \v\ < 1 
leads to the same conclusions. This is one of the main 
results of this paper. 

To quantify the magnitude of discretization and finite- 
size effects we have computed the first two cumulants of 
the topological charge distribution with the estimators 
given in Eqs. (|12p and The numerical results are 

reported in Table [Til The contributions from the poorly 
sampled tail of the distributions, i.e. bins of the sym- 
metrized histogram populated by less than 10 events, 
have been estimated from the large volume expression 
in Eq. (flO)) as suggested in Ref. [2^. Within our sta- 
tistical errors we do not observe a signal for the higher 
cumulants, and correlation functions of an odd number 
of topological charges are always compatible with zero. 

For the data samples that can be directly compared, 
the values of (Q^) are in very good agreement with the 
results in Ref. [17]. We confirm that discretization ef- 
fects on fgCi, where tq is a low-energy reference scale 
well measured in the pure gauge theory [3^, are moder- 
ate (of the order of 10% at our coarser lattice spacing), 
and finite size effects are below 5% at our smaller vol- 
ume. Even though our errors for (Q^) are much smaller 
with respect to those in Ref. , a significative improve- 
ment of the determination of Ci in the continuum limit 
requires further simulations and is left to a future publi- 
cation. 

The second cumulant is best expressed by the adimen- 
sional ratio (Q*)'^°"/(Q^) which has a well defined contin- 
uum and infinite- volume limit. The numerical results for 
the lattices A1-A3 and Bi are reported in Table |TT1 and 
they are plotted in Fig. [5] as a function of (a/rp)^. All 
values are incompatible with 1 , the predicted value from 
Eq. The results from Ai, A2 and A3 agree within 
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(Q > 


{Q } /{Q ) 


Ai 


1.637(13) 


0.60(9) 


0.37(6) 


Aa 


1.566(13) 


0.47(9) 


0.30(6) 


A3 


1.432(12) 


0.43(7) 


0.30(5) 


Bi 


3.09(3) 


0.8(3) 


0.27(10) 


Ci 


5.44(8) 







TABLE II: Results for the first two cumularits and their ratio. 
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FIG. 2: Ratio of the first two cumulants vs the lattice spacing. 



errors. No statistical-significant evidence of discretiza- 
tion effects is tlius observed. The tlieoretical arguments 
given in the third section suggest small finite-size effects 
at these volumes. A direct estimate of these effects on 
{Q'^y°'^ would require more precise data for the lattice 



Ci. The compatibility of the results on lattices Ai and 
Bi, however, is consistent with the theoretical expecta- 
tions. Our best estimate for the ratio of the first two 
cumulants is {Q'^Y°'' / {Q"^) = 0.30(11). The central value 
is taken from the lattice A3, the one with the finer lat- 
tice spacing, and the error is the sum in quadrature of the 
statistical error and of the difference between the central 
values computed on the lattice Ai and Bi. 

Final remarks. — From the previous analysis we con- 
clude that our results disfavour the 9 behaviour of the 
vacuum energy predicted by instanton models. Our best 
estimate of the ratio (Q'*)'=°"/((3^) = 0.30(11) is incom- 
patible with 1, which is the value predicted from Eq. Q. 
This suggests that the quantum fluctuations of the topo- 
logical charge are of quantum non-perturbative nature 
in the ensemble of gauge configurations that dominate 
the path integral. The large Nc expansion does not pro- 
vide a sharp prediction for the value of {Q'^Y°'^ / (Q^). Its 
small value, however, is compatible with being a quantity 
suppressed in the large Nc limit. The value of ((3*)'^°" is 
related via the Witten-Veneziano mechanism to the lead- 
ing anomalous contribution to the rj'-rj' elastic scattering 
amplitude in QCD. 
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